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Abstract 

The pairing of fermionic atoms in a mixture of atomic fermion and boson gases at zero temper- 
ature is investigated. The attractive interaction between fermions, that can be induced by density 
fluctuations of the bosonic background, can give rise to a superfluid phase in the Fermi component 
of the mixture. The atoms of both species are assumed to be in only one internal state, so that the 
pairing of fermions is effective only in odd-/ channels. No assumption about the value of the ratio 
between the Fermi velocity and the sound velocity in the Bose gas is made in the derivation of 
the energy gap equation. The gap equation is solved without any particular ansatz for the pairing 
field or the effective interaction. The p-wme superfluidity is studied in detail. By increasing the 
strength and / or decreasing the range of the effective interaction a transition of the fermion pairing 
regime, from the Bardeen-Cooper-Schrieffer state to a system of tightly bound couples can be 
realized. These composite bosons behave as a weakly-interacting Bose-Einstein condensate. 
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I. INTRODUCTION 



Trapped mixtures of ultracold atomic fermions and bosons in gaseous state offer a con- 
venient testing ground for many-body theories These systems display a rich phase 
diagram. Depending on the strength of the interactions among constituents and on the 
density of the components, several different phases, like a single mixed phase, a coexistence 
of pure and mixed phases, or a collapse of the mixture at densities above some critical value. 



can occur 



Another interesting feature of these systems is that an effective fermion- 
fermion interaction, mediated by the bosonic component, can be induced by the exchange 
of virtual phonons. This fact has been pointed out long ago for dilute solutions of ^He in 



superfluid ^ife [ll| . The induced interaction becomes particularly relevant in ultracold and 
diluted gases when the spins of the fermions are polarized by an external magnetic field since 
in this case the bare interaction between fermions becomes ineffective. This is a consequence 
of the Pauli principle. An interesting property of the induced interaction between fermions 
is that it is attractive, irrespective of the sign of the fermion-boson scattering length, and 
this may lead to the onset of odd-/ superfluid phases in the fermion com pon ent of the mix- 
ture. These fermionic phases would coexist with the bosonic condensate 12|-[l8|. Moreover, 
as pointed out in 1^, results obtained in the study of Bose-Fermi mixtures, could be of 
interest also in the field of dense QCD systems containing bosonic di-quarks and fermionic 
unpaired quarks as effective degrees of freedom. 

The possibility of controlling the strength of the interatomic interactions via Fano- 
Feshbach resonances has given further impulse to the study of these s yste ms. After the 
observation of Fano-Feshbach resonances in Bose-Fermi gas mixtures 2Ch23| , several exper- 



imental investigations of the behaviour of these systems as a function of the inter-species 
interaction strength have been performed [23-27|. On the theoretical side, the study of 
these systems has addressed mainly the problem of stability against collapse, or that of 
phase separation with the possible formation of composite fermionic or bosonic molecules. 
The occurrence of boson-fermion pairing correlations has also been studied ^i^. More- 
over, the authors of Ref. [3] suggest the occurrence of a p-wave superfluid phase in the 
fermionic component of a Bose-Fermi gas mixture when the inter-species interaction be- 
comes sufficiently strong. They further point out the possibility of a transition from long- 
ranged pairing to tightly bound pairs of fermions when the interaction strength increases. 
In Ref. [l81], uniform gases in a spin-polarized state are considered, so that the attractive 
fermion-fermion interaction is provided by the exchange of virtual phonons. 

The evolution of superfluidity from the Bardeen, Cooper, Schriffer (BCS) regime to the 
Bose-Einstein condensation (BEC) limit has been studied for several different systems of 
fermions, including ultracold Fermi gases and nuclear matter. In ultracold Fermi gases, 
the transition occurs when the scattering length changes sign across a Fano-Feshbach reso- 

for p-wave superfluidity), while in nuclear matter, a BCS-BEC 



371 ] (see also [38|, |39 



nance ^ ^ ^ ^ ^ 

transition is expected when the system becomes rarefied 

In this work we want to study the transition from the BCS to the BEC regime for the 
fermionic component of a gaseous mixture of ultracold bosons and spin-polarized fermions. 
Our treatment is at zero temperature, moreover we assume that the system is sufficiently 
diluted so that the bare fermion-fermion interaction can be neglected. An attractive in- 
teraction between fermions arises due to the exchange of virtual phonons and, because of 
the Pauli principle, this interaction is effective only for odd values of the relative angular 
momentum. This induced interaction is derived by extending the results of Ref. [16|] to the 



2 



limit of zero temperature. In 16|] an explicit equation for the energy gap has been derived 
and solved in the limit in which the sound velocity cs in the Bose gas is much smaller than 
the Fermi velocity vp- Here that approximation is abandoned, since densities and coupling 



constants are varied over a wide range of values. In Ref. [18| instead, the opposite limit 
vp has been assumed, but that approximation is also not appropriate for the same rea- 
sons. The obtained equation for the energy gap is solved without introducing any particular 
ansatz ioi the functional form of the pairing field or the effective interaction. Moreover, the 



angle-averaging approximation, used in [16|, for the quasiparticle energies is dropt. This 



permits a better knowledge of symmetry properties of the pairing field. 



II. FORMALISM 



In this paper we rely on the results of Ref. 16[. There, a procedure for determining 
the induced interaction between fermions has been proposed, by using methods based on 
quantum field theory at finite temperature. 

By treating both components of the mixture in mean-field approximation, the resulting 
effective interaction was written as 



Veff{q, Ti - T2) = Vi{q, Ti - T2) + V2{q, n - T2) . 



The first term 



Vi{q,n - T2] 



B 



Di^'\q,n-T2) + Dr{q,r2-n) 



(11), 



+ 2Df^(g,ri-r2) 

is the contribution from states containing one virtual phonon, while the second one 

f - 

(27r)3 J 

+ Di''\k,n-r2)Di,''\\k-ci\,r2-n] 



(1) 



V2{q,Ti - T2] 



/^r(A:,ri-r2)Z^r(|k-q|,r2-ri 



(2) 



is the contribution from states containing two virtual phonons. States with more than two 
virtual phonons have been neglected. 

In the equations above n% is the density of the Bose condensate, A = 2-naBF /thr is the 
boson-fermion effective coupling constant, while a^i? is the fermion-boson scattering length 
and m/j = mBmp / {irtB + nip) is the reduced mass for a boson of mass and a fermion 
of mass mp. The quantities Dq^^ are the components of the imaginary-time Bogoliubov 
propagator (see e.g. sect. 55 of Ref. El for their explicit expressions, note however that 
the definition used here and in [16|| differs by an overall minus sign with respect to that of 

The equation for the pairing field at zero temperature can be obtained by analytic con- 
tinuation to real times (r — )■ it) of Eq. (16) in Ref. 16|. It takes the form 



Afk,k',t-n 



(27r)^ 



j dkirfks K//(ki - k, t - t')^(k + k' - ki - ka) 

xG^''\\i2M,t - t') , (3) 
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where G'(^^)(k2, ki, t — t') is the anomalous propagator for fermions interacting with the 
pairing field A(k, k',t — t') (units h = c = 1 are used). The pairing field is antysimmetric 
under the exchange k k', this is a consequence of the Pauli principle for a couple of 
fermions in the same spin state. 

For vanishing temperature the contribution to the effective fermion-fermion interaction 
from the exchange of two phonons can be neglected [3]. The quantum depletion of the 
boson condensate will also be neglected since, for the purposes of the present work it does 
not play a significant role. Then the density of the Bose condensate coincides with the actual 
density of the Bose gas 

/^B — Anp 
nB = n"j^= . (4) 

7 

Here 7 = AiraBB/^B is the boson-boson effective coupling constant, while ass is the boson- 
boson scattering length. Within a mean-field approach the effect of the interaction with the 
Fermi gas simply amounts to replacing the chemical potential with the effective value 
— ^np, where np is the fermion density [l6| . 

The one-phonon-exchange effective interaction between fermions reads 

Viiq,h - t2) = -X^riB [D'i^\qM - ^2) + D^^^\qM - h) + 2D^^''\q,t^ - t^)] , (5) 

where Dq^^ are now the components of the real-time Bogoliubov propagator; for their explicit 
expressions see e.g. sect. 21 of Ref. 0. In addition the fermions interact with the 



condensate mean field acquiring the energy Xub- This term can be added to the fermionic 
chemical potential, giving an effective potential ii*p = — ^^b to he determined by fixing 
the density of the Fermi gas. 

The solutions of Eq. ([3]) with center of mass momentum P = k + k' = ki + k2 7^ cor- 



respond to the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase j47j. Here, only solutions 
with P = are considered. In this case, the pairing field and the anomalous propagator 
depend on only one momentum, k = — k' and ki = — k2, respectively. In the frequency 
representation Eq. ([3]) reads: 

A(k,u;) = / dkWVi(k - k', w - uj')G^^^\\i:,uj') . (6) 

(zvr)'* J 

In order to simplify calculations the frequency dependence of the pairing field will be ne- 
glected and the static limit A(k, = A(k,u = 0) = A(k) will be taken. This approxima- 
tion amounts to assuming an instantaneous pairing field, neglecting retardation effects for 
the fermion pairing. In the time-dependent representation: 

A(k, t-t') = 6{t - f) J d{t - t') A(k, t-f). (7) 

In this work we are mainly interested in the spatial correlations of paired fermions, and 
this approximation should not change the main results of our calculations, at least at a 
qualitative level. 

When the dependence on uj of the pairing field is neglected, the anomalous Green's 
function is given by the simple expression 
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with ^(fc) = k'^/2mF — fi*p. By replacing this expression into Eq. ([6]) and performing the 
integration over the frequency u' the following equation for the static limit of the pairing 
field is obtained 

""^^^ (2vr)3 J [E{k') + u{q)] Wq)E{k') ) ' 
where e(g) = g^/(2mB), (X'(g) are the excitation energies of the Bose gas calculated 

within the Bogoliubov approximation, u){q) = \J (g^/ (2mB))^ + 2'ynB / {2mB), with q = 
k — k'. The quantities E(k') are the energies of the fermionic quasiparticles, E{k.') = 



^^(/c') + A(k')A*(k'). We notice that the approximation cs ^ vp made in Ref. |18| . 
amounts to neglecting the quasiparticle energies -E'(k'), with respect to the boson energies 
uj(q) in the term in square brackets of Eq. (Q, whereas the opposite approximation of Ref. 
[l6[ is equivalent to assuming E(k') ^ u;(g). 

From the rotational invariance of the coupling between k and k' it ensues that if A(k) 
is a solution of Eq. OH]) also A(7^k), where 71 represents any rotation in the k-space, is a 
solution. This implies a very high degeneracy for the ground-state energy of the superfluid 
phase. We will discuss this point below. 

The field A(k) can be expanded in spherical harmonics, 



4-71 

with only odd values of / contributing because of the antisymmetric relative-motion wave- 
function. In Eq. (Q components of A(k) with different / are coupled in general. We have 
estimated the relative weight of the components with / > 1 using an angle-average ap- 
proximation for the quasiparticle energy [l6|. In this approximation the /-components are 
decoupled, and in all the cases studied in the present work those with / > 1 turn out to be 
very small with respect to the / = 1 component. Hence, we restrict our investigation to the 
I = 1 component alone. Explicit calculations show that this approximation is fairly reliable 
when we limit ourselves to evaluate average quantities, that is, quantities integrated over 
the relative momentum of the couple of fermions. 

Taking into account only the I = 1 component, the pairing field can be put as A(k) = 
^i{k) ■ k, with k = k/k. From Eq. (jO]) one can easily check that the three-component 
quantity ^i{k) behaves as a vector. Then, the superfluid state can be identified by this 
vector. One also can see that if Ai(fc) represents a solution of the gap equation also 
TZA.i{k), where TZ denotes an arbitrary rotation in the orbital coordinates, is a solution. 
This means that only the magnitude Ai(/c) can be fixed by the gap equation, whereas the 
direction of Ai(/c) is completely arbitrary. Finally we remark that the energy gap vanishes 
when k hes in the plane perpendicular to Ai(/c). 

The expression for the ground-state energy of the Fermi gas in the superfluid phase can 
be derived from Eq. (15) of Ref. [l6| for the effective action of the pairing field, by adding 
the contribution of a noninteracting Fermi gas with the effective chemical potential fip. By 
exploiting the equation for the pairing field (Eq. (16) of Ref. [l6|) and taking the limit of 
vanishing temperature the energy per fermion in the superfluid phase is given by 
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This expression coincides with the usual expression of the BCS theory (see, e. g., Ref. |46|). 
apart from a factor 1/2 due to the absence of degeneracy for the Fermi gas in the present 
case. 

In order to determine exphcity the field Ai ■ k together with the effective chemical po- 
tential fip the equation fixing the fermion density has to be added 



Hp 



fj^u,_m\. (11) 

J (27r)3 2 I E(k) J ^ ^ 



Both Eqs. f IlOp and fllip are invariant under the transformations Ai(A;) — > TZA.i{k). 
Then, the ground-state energy and the chemical potential of the superfiuid phase are deter- 
mined by the magnitude Ai(A;) alone. This implies that the energy of the superfiuid ground 
state is infinitely degenerate. Fermions can test states of the degeneracy subspace through 
rotations of the vector Ai(A;). We notice that in the present case the spins of fermions do not 
play any role since they are frozen along a fixed direction and there is no coupling between 
spin and orbital degrees of freedom. 

In order to calculate its magnitude we choose a particular direction for the vector Ai(/c), 
say Ai(fc) = (0, 0, Ai(fc)). Moreover, we observe that the scaled quantity, Ai(A;) = Ai{k)/eF 
where ep = k\/2mF = {^'n'^^n'l'^ jimp is the Fermi energy, depends only on three dimen- 
sionless quantities: 

h = {\/'^){\np I tp), c = ^Bkp, d = mB/mp. 

The parameter h is determined by the mean field of fermions acting on bosons and by the ratio 
between the fermion-boson and boson-boson coupling constants, in practice it represents 
the strength of the effective fermion-fermion interaction. The parameter c is the coherence 
length of the Bose condensate = l/\/2mBTraB in units of l/kp, hence it represents the 
ratio of the range of the effective interaction between fermions to the average interparticle 
spacing of the Fermi gas. In terms of these dimensionless parameters the equation for Ai(fc) 
reads 

Ai(fc) = j^b I dn^cos{e^) I -^cos(ek')^^ \ . , ^'^^'^ 



(12) 

with -E'(k') = _E'(k')/ei?, ^{k) = ^{k)/ep, and the momenta are expressed in units of the 
Fermi momentum: k = k/kp and k' = k'/kp. 



III. RESULTS 

In this section the properties of the fermionic component of the Bose-Fermi mixture as 
functions of the two dimensionless parameters b and c are examined and discussed. 

The value of the parameter d = mB/mp is fixed by specifying the components of the par- 
ticular system considered: a mixture of ^^Rb (bosons) and '^°K (fermions) atoms. Once the 
ratio of the masses is fixed, the ratio between the Fermi and phonon velocities is determined 
by the value of the parameter c alone, since vp/cs = \/2cd. 

Physical quantities of interest will be studied both as a function of the parameter c for 
a fixed value of 6 (6 = 5) and as a function of b for a fixed value of c (c = 2). The first 
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FIG. 1: The maximum value of the p-wave energy gap in units of the Fermi energy as a function 
of 1 /c with 6 = 5 (circles) , and as function of b with c = 2 (squares) . 



case can be implemented by varying for instance the boson density alone, whereas for the 
latter the boson-fermion coupling constant can be varied while the other parameters are 
kept constant. The range of values chosen for the parameters b and c includes domains 
where the properties of the Fermi gas change drastically. The peculiar behaviour of various 
physical quantities generally suggests the occurrence in the Fermi gas of a transition from 
a superfiuid long-ranged phase (BCS phase) to a phase of tightly bound pairs of fermions. 
These composite bosons behave like a weakly interacting Bose condensate (EEC phase). For 
6 = 5 and c = 2 the obtained energy gap is very small, implying that the fermions are well 
inside the BCS phase. 

It should be remarked that in the considered range of values for c the ratio of the Fermi 
velocity to the phonon velocity lies in the range 1.2 < vp/cs ^ 9, so the approximations 



cs ^ vp and cs <C vp, used in Ref. [18| and in Ref. [16[ respectively for deriving the gap 
equation, are not valid in the region close to the BCS-BEC transition. 

In Fig. [T]the maximum value of the magnitude of A.i{k), Ai^ax, in units of ep, is shown 
both as a function of c and as a function of b. The value of Aimax increases with increasing 
b or decreasing c, i.e. when the effectiveness of the fermion-fermion interaction induced by 
the exchange of phonons increases. A similar behaviour is shown by the value of k at the 
maximum of Ai(A;), i.e. the peak of Ai(fc) moves from the Fermi surface towards higher 
values of k within the range ~ {kp,4kp). We can also see that Aimax grows quickly with 
1/c for 1/c ^ 1.5, while it shows a smoother behaviour as a function of b instead. Finally, 
it should be remarked that, as an odd-/ component, Ai(A;) vanishes for A; — )■ 0, contrary to 
the case of s-pairing. 

Figure |2] shows the behaviour of the effective chemical potential fi*p and of the energy 
per particle Ep for the Fermi component of the mixture as a function of 6 or 1/c. For low 
values of b and/or high values of c the relations fi*p ~ ep and Ep ~ 3/5 ep hold according to 
the weak coupling BCS theory. When the effectiveness of the interaction between fermions 
increases, fi*p and Ep start to decrease becoming eventually negative. In addition, their 
values become similar when the fermion-fermion coupling increases. These peculiar features 
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FIG. 2: Effective chemical potential (filled red symbols) and energy per particle (empty blue 
symbols) of the fermionic gas, in units of the Fermi energy, as a function of 1/c with 6 = 5 (circles), 
and as function of b with c = 2 (squares). 



indicate the formation of bound pairs of fermions behaving as a condensate Bose gas (BEC 
phase) [13]. In particular, the difference ^y — Ep-, which is related to the interaction between 
the composite bosons, approaches the values 0.02 for c = 0.35 (with 6 = 5) and 0.06 Ci?' 
for 6 = 60 (with c = 2). 

A further remark is required. Equation ([H]) for the pairing field corresponds to the saddle- 
point condition for an effective action [l6|. Thus, our calculations are essentially performed 
within the framework of a mean-field theory. The extension of our approach to situations 
where the energy gap can be larger than the Fermi energy might seem inappropriate. How- 
ever, analogous studies 3^, 45 1 on the superfluidity in Fermi gases have shown that, provided 



the temperature is much lower than the critical one, fluctuations of the pairing field about 
the saddle-point value do not play an important role in determining the energy gap. Hence, 
we confide that our results can be correct, at least qualitatively. 
We turn now to the momentum distribution of the fermions 

The anisotropy of the pairing field changes the features of n(k) significantly with respect to 
the s-wave superfluidity case. Choosing the direction of Ai(A;) along the z-axis the energy 
gap vanishes for k lying on the plane = 0. Accordingly n{kx, ky, 0) is given by the step 
function 9{—C,{k)) on the BCS side { fi*p > ) and vanishes on the BEC side { fip < ). 
In the BCS phase the fermion distribution is concentrated across the plane k^ = 0. In the 
BEC phase instead, fermions fill two domains, which are symmetrical with respect to the 
plane k^ = 0. The angle-averaged distribution still shows a critical behaviour when the 
fermionic component of the mixture crosses the borders between the BCS side and the BEC 
side. This effect can be observed in Figs. [3] and HI where the momentum distribution is 
displayed for two sets of values of b and c, which include both the sides of the transition. 
On the BCS side the step of the momentum distribution becomes less and less pronounced 
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FIG. 3: Main figure: angle-averaged momentum distribution of tlie fermions, n{k), on the BEG 
side for two different values of c with fixed b = 5: c = 0.45, solid line, and c = 0.4, dashed line. 
Inset: same as in main figure but with c = 0.6, solid line, and c = 0.48, dashed line (BGS side). 
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FIG. 4: Main figure: angle-averaged momentum distribution of the fermions, n{k), on the BEG 
side for two different values of b with fixed c = 2: b = 35, solid line, and b = 50, dashed line. Inset: 
same as in main figure but with b = 25, solid line, and b = 30, dashed line (BGS side). 



as becomes smaller. On the BEG side the distribution shows the typical behaviour of 



an odd-/ pairing: it vanishes together with Ai{k) at fc ^ |l8l . |39| . Moreover its peak 
becomes less pronounced and moves toward higher value of k/kp, when the effectiveness of 
the induced interaction between fermions increases. 

Here, we do not discuss the effects of the anisotropy of the pairing field on the quasiparticle 
spectrum, i.e. a gapless to gapped quantum phase transition, when the effective chemical 
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FIG. 5: Main figure: radial wave-function with / = 1 of a pair of fermions on the BEC side as 
a function of the ratio v/tq (rg is the average spacing of the fermions) for two different values of 
c: c = 0.45 black line and c = 0.4 green line, with 5 = 5 in both cases. The wave-functions are 
normalized according to Eq. (jl3p . Inset: tail of the pair wave-functions (solid lines) and the wave 
functions of a particle in a square-well potential (dashed lines), see text. 



potential vanishes /x^ = 0, a topic which has been widely discussed in Refs. 18|, |38|, |39 



We only observe that in the BCS phase the quasiparticle energy vanishes when k is on the 
Fermi surface, ^{k) = 0, and perpendicular to Ai(A;). 

The transition towards a gas of bound pairs of fermions can be better appreciated by 
looking at the spatial structure of the pairs. The pair wave function in the center of mass 
frame is obtained from the Fourier transform of the anomalous density 



/Ci(k) 



Ai(A;)-k 
2E{k) ■ 



With our choice for the direction of the vector Ai(fc) the pair wave function can be written 

as 

with only odd values of /. 

The wave function 0(r) can contain several partial waves. However, explicit calculations 
show that the norm of the I = 1 component exhausts the norm of 0(r) within a few percent. 
Hence, we can neglect the components with I > 1. 

It is convenient to normalize the radial wave-function (f)ifi{r/ro) according to the condi- 
tion 



d 



(13) 



where tq is the average spacing of the fermions. On the BCS side the radial wave-function 
shows the usual damped oscillatory behaviour, where the first peak becomes larger than the 
secondary oscillations as fi*p approaches zero. On the BEC side the oscillations disappear 
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FIG. 6: Same as in Fig. [5] but for 6 = 35 black lines and 6 = 50 green lines, with c = 2 in both the 
cases. 

instead, and the wave function acquires a shape similar to that of a bound state with 1 = 1. 
In Figs. 13 and E] the radial wave-function is displayed as a function of r/rg for two different 
values of the parameters b and c. We can see that when increasing b and/or decreasing 
c the wave function is squeezed within a narrower domain about the origin. This simply 
means that increasing the strength and/or decreasing the range of the induced interaction 
the binding of a pair of fermions becomes tighter. 

The insets of Figs. [5] and [6] show a comparison of the asymptotic behaviour of the radial 
functions for a pair of fermions with that of the radial wave functions for bound states 
with / = 1 of a particle in a square well potential. For the binding energy and the mass 
of the bound particle we take the energy and the reduced mass of the pairs of fermions, 
BE = 2\Ef\ and mji = m/2. The point where two wave functions with the same binding 
energy have been joined, is far enough from the position of the peak of the pair wave-function 
and is assumed to be outside the potential well. Once more, this comparison suggests the 
occurrence of a transition of the Fermi component of the mixture toward a Bose gas of 
tightly bound dimers, when the effectiveness of the induced interaction increases. The root 
mean square radius of the dimers approaches the values ~ 0.2ro for c = 0.35 (with 6 = 5) 
and ~ O.Stq for 6 = 60 (with c = 2). This implies that the density of the couples of fermions 
practically is half that of the Fermi gas. 

The calculated values of the difference fi*p — Ep allow us to give an estimate of the ratio 
between the effective scattering length and the interparticle spacing for the bound pairs 
of fermions. In the ladder approximation for the self-energy of a dilute Bose gas [46| with 
repulsive interaction, the relation between the chemical potential and the energy per particle 
is given by 

Q4_ln^y (14) 

where n^, ad and respectively represent the density, scattering length and mass of the 
constituents of the Bose gas, with Ud = np/^, rrid = 2mF and /i — E/N = 2(yU* — Ep) for 



E 2Txndad 



N rrid 
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the composite bosons. Dividing Eq. ( !T^ by the Fermi energy one obtains 




t * T7 \^ 1 1 1/3 

With the calculated values for the l.h.s. of the equation above the product a^nd approaches 
the values 1.7 10^^ for c = 0.35 (with 6 = 5) and 1.5 10"^ for 6 = 60 (with c = 2). These 
values are sufficiently low to allow us to consider the system of composite bosons as a weakly 
interacting Bose-Einstein condensate. 

The repulsive nature of the force between composite bosons of the BEC phase is consistent 
with the sign of the difference between the chemical potential and the energy per particle 
for fermions in a superfluid phase. In fact Eq. (fTOl) clearly shows that anyhow fi*p — Ep 
is positive. On the other hand a repulsive interaction is required for the stability of the 
BEC phase against collapse. Our result to some extent is in agreement with the studies of 
Refs. Blii on the few-body properties of p-wave molecules of fermionic atoms. However, 



the calculations of Refs. do not rule out the possibility of an attractive interaction 

between composite particles. 



IV. SUMMARY AND CONCLUSIONS 



The transition to the BEC regime of a superfluid Fermi gas in a Bose-Fermi mixture 
can be obtained by tuning the strength of the interatomic interactions via Fano-Feshbach 
resonances and/or by varying the densities of the two gases. In this paper the relevant 
quantities for the transition have been expressed as functions of two mutually independent 
dimensionless parameters h and c. These two parameters are related to the strength and 
range of the effective interaction between fermions induced by the exchange of one virtual 
phonon. When the strength of the induced interaction increases and/or its range decreases 
the fermion pairing evolves from a situation of long-ranged correlations to the onset of 
tightly bound states. The occurrence of bound states of couples of fermions can be further 
evinced by the asymptotic behaviour in coordinate space of the anomalous density. This 
quantity can be interpreted as the relative-motion wave function of the couples of fermions 
and its asymptotic behaviour approaches that of a particle bound in a spherical potential 
well. Moreover, when the effectiveness of the induced interaction increases, the chemical 
potential and the energy per particle of the fermionic component of the mixture approach 
each other. This behaviour suggests that the pairs of fermions behave as a Bose-Einstein 
condensate with a weak repulsive interaction between the composite bosons. Finally, the 
root-mean square radius of the bound pairs is a small fraction of the average interparticle 
spacing of the fermion gas, so that the gas of composite bosons practically has the same 
degree of diluteness as the Fermi component of the mixture. 
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